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ABSTRACT  . 

An  investigation  was  conducted  into  the  possibility 
of  visually  detecting  a  dewetted  zone  in  a  solid  rocket 
propellant  through  the  use  of  holographic  interferometry , 
Techniques  used  were  simple  strain,  residual  strain  and 
time-average  holography.   As  a  basis  for  comparison  with 
experimental  data  a  finite  element  program  was  used  to 
model  the  test  samples. 
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I.   INTRODUCTION 

A.  DEWETTING 

Dewetting  is  the  phenomenon  that  makes  filled  materials 
(rocket  propellants)  differ  from  classical  viscoelastic 
(unfilled)  materials.   Propellant  consists  of  a  composite 
of  solid  oxidizer  particles  imbedded  in  and  bonded  to  a 
fuel  matrix.   After  an  initial  linear  region,  nonlinear 
stress-strain  behavior  under  loading  results  from  the  pro- 
gressive failure  of  these  bonds  at  the  interface  between 
the  fuel  matrix  and  the  oxidizer  particles.   Voids  or 
vacuoles  are  thus  created  throughout  the  materials  and  their 
volume  increases  with  further  loading,  causing  an  overall 
volume  increase  in  the  material.   This  phenomenon  is  called 
dewetting. 

B.  DEWETTING  THEORY 

Lindsey  [1]  has  postulated  that  a  dewetting  criterion 
exists  and  that  it  is  expressible  as 

J2  +  (Xd]^)2  =  constant 
where  J2  is  the  second  deviatoric  stress  invariant,  I   the 
first  stress  invariant  and  0(  is  a  small  constant.   This 
expression  represents  a  surface  in  principle  stress  space 
that  is  axisymmetric  about  the  hydrostatic  axis  {O^-Op-O^). 
Figure  1  is  a  sketch  of  the  general  shape  of  this  surface. 
A  state  of  stress  located  outside  of  or  on  this  surface 
will  be  dewetted. 


C.   EVALUATION  OF  THEORY 

The  dewetting  created  by  excess  strain  is  not  visible 
to  the  naked  eye.   However,  changes  in  the  dewetted  material 
properties  will  result  in  a  displacement  field  that  is  dif- 
ferent from  that  of  a  material  which  has  similar  properties 
but  does  not  dewet.   To  test  the  postulate  of  Lindsey, 
holographic  interferometry  was  selected  to  detect  these  dis- 
placement changes  because  (1)  it  is  sensitive  to  very  small 
displacements  and  (2)  holography  models  need  not  be  trans- 
parent and  birefringent  as  in  the  case  of  photoelasticity. 

Attempts  to  identify  a  dewetted  area  through  the  visual 
interpretation  of  the  interf erometric  fringe  patterns  in- 
clude the  techniques  of  small  strains,  residual  strains  and 
vibrational  time  average. 

A  finite  element  program  was  used  as  a  basis  of  compar- 
ison between  the  displacements  obtained  holographically  from 
both  dewetting  and  nondewetting  materials  of  similar  shape 
undergoing  small  strains. 


II.   EXPERIMENTAL  APPARATUS 

A.  THE  HOLOGRAPHIC  SETUP 

The  holographic  setup  was  of  a  conventional  nature. 
An  airtable  was  used  consisting  of  a  two-inch  thick  piece 
of  plywood,  four  by  eight  feet  in  dimensions,  supported  by 
four  inner  tubes  resting  on  a  heavy  metal  table.   Very  good 
vibration  insulation  was  obtained. 

A  Spectra-Physics  15  milliwatt  continuous  wave  helium- 
neon  laser  beam  was  expanded  through  a  spatial  filter  and 
directed  onto  the  sample.   A  front  surface  mirror  was 
mounted  beside  the  sample  to  provide  a  reference  beam.   The 
path  length  from  source  to  sample  to  holographic  plate  was 
approximately  the  same  as  the  path  length  from  source  to 
reference  mirror  to  plate.   The  holographic  plate  was  mount- 
ed directly  in  front  of  the  sample.   Figure  2  is  a  schematic 
and  photograph  of  the  setup. 

The  film  used  was  Agfa-Gevaert  8E75  Scientia  four-by- 
five  inch  plates.   The  developing  procedure  was  five  minutes 
in  Kodak  D-19  developer,  30  seconds  in  a  stop  bath  and  five 
minutes  in  a  fixer  followed  by  a  30-minute  water  bath. 

B.  THE  TENSILE  TESTING  MACHINE 

The  use  of  a  conventional  tensile  testing  machine  was 
rejected  because  it  has  only  one  moving  crosshead.   Strain- 
ing a  sample  from  one  end  results  in  the  center  of  the 
sample  moving  about  half  the  distance  of  the  crosshead. 


This  movement  is  undesirable  in  holographic  work  where  even 
pure  rigid  body  movement  causes  fringes  in  double  exposure 
interf erograms . 

To  minimize  the  translation  of  the  center  of  the  sample, 
a  tensile  testing  machine  was  built  that  had  two  moving 
crossheads.   Turning  the  handwheel  resulted  in  the  cross- 
heads  moving  in  opposite  directions.   One  complete  turn  of 
the  handwheel  resulted  in  each  crosshead  moving  0.005  inch. 
Figure  3  is  a  photograph  of  the  machine. 

The  handwheel  had  notches  in  it  every  five  degrees  which 
engaged  a  springwire  latch.   This  allowed  the  handwheel  to 
be  turned  in  multiples  of  five  degrees  in  the  dark  for 
double  exposure  work. 

C.   THE  TEST  SAMPLES 

The  test  samples  were  six-by-six  inch  squares,  one- 
quarter  inch  thick,  with  a  one-inch  diameter  hole  in  the 
center  to  act  as  a  stress  riser.   A  one-inch  piece  of  wood 
was  bonded  to  two  sides  of  the  samples  to  provide  a  place 
for  mounting  in  the  tensile  testing  machine.   This  type  of 
geometry  was  selected  for  its  ease  of  construction  and  its 
two-dimensional  plane  stress  properties. 

The  sample  material  was  obtained  from  United  Technology 
Corporation,  Sunnyvale,  California.   The  unfilled  material 
was  rocket  motor  case  liner  of  type  C-296-182-7  and  the 
filled  material  was  rocket  propellant  of  type  UTP-3001. 
Figure  4  is  a  drawing  of  the  sample  geometry. 


III.   DISPLACEMENT  FIELD  DETERMINATION 

A.   COMPUTER  GENERATED  DISPLACEMENTS 

1 .  The  Finite  Element  Program 

As  a  basis  for  comparison  with  the  experimental  data 
a  finite  element  program  called  PSELST  (Plane  Stress  Elastic 
Analysis  Using  Linear  Strain  Triangles)  was  used  to  model  the 
test  samples.   This  program  was  able  to  handle  plane  stress 
problems  with  up  to  six  different  types  of  material  in  its 
composition.   This  is  a  necessary  feature  in  order  to  handle 
the  material  property  changes  caused  by  dewetting.   The  out- 
puts of  the  program  were  the  stresses  and  displacements  at 
each  node  in  the  model.   Reference  2  describes  PSELST  in 
great  detail. 

2.  The  Finite  Element  Grid 

The  finite  element  grid  is  shown  in  Figures  5  and  6. 
Only  the  upper  right-hand  quadrant  of  the  sample  was  modeled 
due  to  the  symmetry  of  the  sample  and  of  the  loading.   The 
smaller  elements  were  placed  around  the  hole  where  the 
stress  gradients  would  be  the  highest.   The  origin  of  the 
grid  coordinate  system  was  placed  at  the  center  of  the  hole 
with  the  X-axis  as  the  horizontal  lower  boundary  and  the  Y- 
axis  as  the  vertical  left  hand  boundary. 

3.  The  Boundary  Conditions 

The  boundary  conditions  on  the  problem  were  straight- 
forward.  Nodal  points  on  the  X-axis  were  free  to  slide  in 
the  plus  or  minus  X-direction  only  and  the  nodal  points  on 


the  Y-axis  were  free  to  move  in  the  plus  or  minus  Y- 
direction  only.   Nodal  points  along  the  upper  boundary 
were  all  given  a  uniform  positive  Y  displacement  with  no 
X  displacement  to  simulate  a  crosshead  moving  up  and  strain- 
ing the  sample.   The  right  edge  nodal  points  were  a  free 
boundary. 

4.   The  Unfilled  Material  Solution 

The  material  properties  for  the  unfilled  sample  were 
taken  to  be  470  psi  for  the  elastic  modulus  and  0.485  for 
Poisson's  ratio,  which  are  typical  values  for  this  type  of 
material.   Experiments  have  shown  that  this  material  is 
linear  to  about  fifteen  percent  strain.   For  this  reason  it 
was  assumed  for  the  computer  solution  that  the  modulus  of 
elasticity  and  Poisson's  ratio  were  constants. 

It  can  be  shown  that  for  a  displacement  boundary 
condition  problem,  the  value  of  the  elastic  modulus  does 
not  affect  the  resultant  displacement  field  if  Poisson's 
ratio  is  held  constant.   As  a  sensitivity  test,  two  computer 
displacement  solutions  were  obtained,  one  with  a  modulus  of 
470  psi  and  the  other  with  a  modulus  of  1090  psi.   Both  had 
a  Poisson's  ratio  of  0.485.   The  resultant  displacement 
fields  agreed  to  within  0.0003  wave  lengths  of  light  for  a 
boundary  extension  of  18.906  wave  lengths. 

A  second  sensitivity  test  was  conducted  with  the 
elastic  modulus  held  constant  at  470  psi  and  Poisson's  ratio 
set  at  0.485  for  the  first  run  and  0.495  for  the  second. 
The  change  in  the  displacement  field  v/as  small  and  probably 
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within  the  error  range  of  the  holographic  technique.   For 
example,  for  a  boundary  extension  of  24.136  wave  lengths 
of  light,  the  difference  between  the  predicted  X  displace- 
ments was  0.0579  wave  lengths  at  grid  point  33. 

The  results  of  the  sensitivity  tests  indicated  that 
the  computer  displacements  for  the  unfilled  material  were 
unaffected  by  an  elastic  modulus  change  of  every  element  if 
Poisson's  ratio  was  held  constant  and  that  the  solution  was 
relatively  insensitive  to  a  Poisson's  ratio  change  of  every 
element  if  the  elastic  modulus  was  held  constant. 
5.   The  Filled  Material  Solution 

The  modulus  of  elasticity  for  the  filled  material 
was  obtained  from  a  uniaxial  tensile  test.   (Figure  7). 
The  stress-strain  curve  was  linearized  by  dividing  it  into 
four  segments  and  replacing  each  segment  with  a  straight 
line.   Poisson's  ratio  for  the  filled  material  was  estimated 
from  data  obtained  by  Lindsey  and  Wood  [1] .   Table  I  is  a 
summary  of  the  data.   The  departure  of  the  stress-strain 
curve  from  linearity  is  due  to  the  onset  of  dewetting  and 
the  stress  at  that  point,  00   ,  is  called  the  yield  or  de- 
wetting  stress. 

The  first  segment  of  the  curve  starting  at  zero 
stress  was  considered  the  nondewetted  portion  and  was 
designated  material  of  type  one.   Type  two  material,  which 
was  partially  dewetted,  was  represented  by  the  second  seg- 
ment of  the  curve.   Type  three  was  dewetted  to  a  greater 
degree  and  type  four  was  considered  to  be  completely  dewetted 
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The  upper  end  of  the  fourth  segment  was  the  stress,  0mnx      , 
at  which  probable  failure  of  the  material  would  occur. 

To  relate  the  two-dimensional  stress  of  the  computer 
solution  to  the  one-dimensional  stress-strain  curve,  an 
equivalent  stress  0       was  used.   Equivalent  stress  [3]  is 
defined  as 


0"e  =-\/3J2 

j  o  =  if(crx-  cp+(cx-cr2)2+  (tfy-  cz)2 


+6(rxy+rxVry2z 

For  plane  stress 

tfz  =  rx2  =  ryz  =  o 

Thus 

Oe =  (o;2  +cry2  -crx0y  +3rx2y)" 

This  equivalent  stress  was  used  because  the  value  of  OC   has 
been  found  to  be  very  small  for  this  type  of  material.   For 
typical  values  of  CX  and  the  stress  field  of  the  specimen, 
the  dewetting  criterion  and  the  equivalent  stress  agree  to 
within  one  percent. 

For  the  computer  solution  the  upper  edge  was  extend- 
ed 0.01  inches  at  a  time.   Although  the  overall  properties 
are  nonlinear,  each  0.01  inch  step  was  considered  small 
enough  to  be  linear.   After  each  extension,  the  output 
stresses  0X,  Oy  and  T^y  at  each  node  were  used  to  determine 
the  equivalent  stress  at  that  node.   If  the  equivalent 
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stress  was  less  than  00    ,  the  node  was  considered  to  still 
be  of  material  type  one.   When  the  equivalent  stress  fell 
within  the  range  of  segment  two,  the  node  was  considered  to 
be  material  of  type  two  and  so  on.   When  the  equivalent 
stress  exceeded  Opnax  /  failure  was  probable  at  that  node. 
Each  quadrilateral  finite  element  had  nine  nodes  associated 
with  it,  but  only  four  were  nodal  points  in  the  grid.   The 
program,  internal  to  itself,  added  a  node  to  the  midpoint 
of  each  side  and  at  the  intersection  of  the  element's  diag- 
onals.  It  was  decided  that  whenever  any  two  of  the  element's 
nine  nodes,  six  in  the  case  of  a  triangular  element,  had 
changed  to  the  next  higher  material  type,  then  the  element 
itself  would  be  changed  to  that  next  higher  material  type. 
Once  an  element  or  node  reached  material  type  four  it  re- 
mained type  four. 

Originally  the  criterion  was  that  five  of  nine  nodes 
had  to  change  to  the  next  highest  material  type  before  an 
element  was  changed  to  that  type.   This  resulted  in  a  rapid 
stress  buildup  at  node  one  and  failure  was  predicted  after 
only  0.11  inch  of  upper  edge  extension.   Thus  a  less  strin- 
gent element  changing  criterion  was  adopted  in  an  attempt 
to  relieve  the  stress  buildup  at  node  one,  because  it  was 
known  from  experiments  that  the  actual  samples  did  not  fail 
until  0.225  inches  of  extension. 

In  addition  to  checking  for  material  changes,  the 

incremental  nodal  displacements  were  added  to  the  nodal 

coordinates  to  obtain  the  new  coordinates  for  the  next 

extension. 
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The  incremental  procedure  used  to  extend  the  comput- 
er solution  to  0.20  inches  may  be  summarized  as  the 
following  iterative  routine: 

a.  Extend  the  upper  boundary  0.01  inch. 

b.  Add  the  incremental  stresses  to  a  running  total 
of  stresses  from  the  previous  extensions. 

c.  Compute  the  nodal  equivalent  stresses  based  on 
the  summation  of  stresses  and  assign  each  node 
its  proper  material  type  based  on  the  linearized 
stress-strain  curve. 

d.  Check  each  element  to  see  if  any  two  of  its 
nodes  have  changed  to  the  next  higher  material 
type.   If  yes,  change  that  element's  properties 
to  the  next  higher  type. 

e.  Add  the  incremental  displacements  to  the  nodal 
coordinates  to  form  the  new  nodal  coordinates 
for  the  next  incremental  extension. 

f.  Return  to  step  a. 

Using  the  above  technique,  failure  was  predicted  to 
occur  at  node  one  after  only  0.12  inches  of  upper  edge  ex- 
tension.  It  was  known  that  the  samples  were  extended  to 
0.225  inches  before  actual  failure,  so  the  computer  predic- 
tion of  failure  was  ignored  and  the  solution  continued  on  to 
0.20  inches  of  extension. 

If  the  dewetting  criterion  is  applied  to  the 
classical  solution  for  the  stress  field  around  a  hole  in 
an  infinite  sheet,  the  outline  of  the  dewetted  zone  will  be 
as  shown  in  Figure  8.   For  comparison,  the  computer  solution 
for  the  dewetted  zone  in  the  filled  material  using  equivalent 
stress  is  shown  for  0.13  inches  of  extension  in  Figure  9. 
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6 .   Nondimensional  Displacements 

The  X  and  Y  displacements  obtained  from  PSELST  were 
due  to  a  0.01  inch  extension  of  the  upper  edge.   However, 
the  unfilled  material  was  considered  linear  so  that  dis- 
placements for  extensions  other  than  0.01  were  obtained  by 
scaling. 

Predicting  displacements  for  the  filled  material  was 
more  restrictive.   The  displacements  had  to  be  scaled  from 
the  X  and  Y  displacements  of  the  twentieth  step  of  the 
incremental  routine. 

All  displacements  were  nondimensionalized  by  divid- 
ing by  the  wave  length  of  light  from  the  Helium-Neon  Laser. 

B.   DISPLACEMENTS  BY  HOLOGRAPHIC  MEANS 
1.   Holographic  Displacement  Theory 

The  method  presented  here  was  developed  by  S.K.  Dhir 
and  J. P.  Sikora  [4].   If  a  light-scattering  object  undergoes 
an  elastic  deformation  between  two  holographic  exposures 
constituting  the  interf erogram,  the  fringes  will  be  related 
to  the  displacement  field.   Let  the  point  P(X,Y,Z)  be  dis- 
placed to  P1 (X+u, Y+v, Z+w)  where  u,  v  and  w  are  the  cartesian 
components  of  the  displacement  vector  (Figure  10) .   Points 
H(Xk,Yk,Zk)  and  S(Xo,Yo,Zo)  are  the  observation  points  on 
the  hologram  plate  and  the  light  source  respectively.   The 
total  distance  traveled  by  the  light  will  be  (SP+PH)  before 
movement  and  (SP'+P'H)  after  movement.   The  difference  in 
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path  lengths  can  be  expressed  as  some  multiple  of  the  wave 
length  of  the  laser  light,  L, 

nL  =  (SP'+P'H)  -  (SP+PH)  (1) 

where  n  is  any  number.   If  n  is  a  whole  number  then  the 
object  point  will  appear  bright  because  the  two  waves  will 
be  in  phase.   If  n  is  some  odd  multiple  of  1/2  then  the 
waves  will  arrive  out  of  phase  and  the  object  point  will 
appear  dark. 

SP1  =  [  (X-Xo+u)  2+(Y-Yo+v)  2+(Z-Zo+w)  2]2" 

SP   =  [  (X-Xo)2+(Y-Yo)  2+(Z-Zo)  2]2"  (2) 

P'H  =  [  (X-Xk+u)  2+(Y-Yk+v)2+(Z-Zk+w)2]  *" 
PH   =  [  (X-Xk)2+(Y-Yk) 2+(Z-Zk)2] 
Rearrange  equation  (1) 

nL  =  (SP'-SP)  +  (P'H-PII)  (3) 

Substitute  (2)  into  (3)  and  the  first  term  becomes 

SP'  =  SP   1  +  _1_   /  2[(X-Xo)u  +  (Y-Yo)v  +  (Z-Zo)w] 

1      (SP)2^        2     2     2\l4 

+  uz  +  vz  +  wz  \ 

Neglecting  small  quantities  in  the  expansion  of  SP1  it  can 

be  shown  that 

(SP'-SP)  =  _l|  (X-Xo)u  +  (Y-Yo)v  +  (Z-Zo)wl        (4) 

SP  L  J 

The  distance  (P'H  -  PH)  can  be  found  by  replacing 

the  subscript  o  with  k  and  1/SP  with  1/PH. 

Now 

nkL  =  A  u  +  B  v  +  C^w  (5) 
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Where 

Ak  =  x"Xo  +  x"xk 
SP      PH 

Bk  =  Y~Yo  +  Y"Yk  (6) 

SP      PH 

ck  =  z~Zo  +  z~zk 
SP      PH 

The  subscript  k  refers  to  a  single  point  of  obser- 
vation on  the  holographic  plate  (Figure  11)  and  n^  is  the 
respective  fringe  number  on  the  point  being  observed. 

The  same  object  point  can  be  observed  through  three 
'   different  points  on  the  plate.   The  result  is  three  equations 
in  three  unknowns  which  can  be  solved  for  u,  v,  and  w. 

If  it  is  not  possible  to  determine  the  location  of 
the  zero  order  fringe,  Equation  (5)  is  modified  to 

(lVnk+l)L  =  (Ak-Ak+l)u  +  (Bk-\+1)v  +  (ck-ck+l>w   (7) 
where  (nk~nk+l^  "*"S  t*ie  num^Gr  °f  fringes  that  shift  through 
an  object  point  going  from  position  k  to  k+1  on  the  plate. 
When  this  subtraction  is  carried  out,  all  the  terms  involving 
1/SP  drop  out.   Thus  when  using  these  equations  with  an  un- 
known zero  order  fringe  the  location  of  the  source  is  unimpor- 
tant.  This  eliminates  errors  in  determining  the  physical 
coordinates  of  the  source.   This  advantage  is  offset  by  the 
necessity  of  combining  two  observations  of  the  same  point  to 
determine  one  equation.   The  sign  assigned  to  a  fringe  shift 
is  not  important  as  long  as  the  system  used  is  consistent. 

In  theory  it  should  take  only  three  observations  to 
determine  the  three  unknowns.   However,  the  value  of  the 


17 


determinant  of  these  equations  is  very  small  due  to  the 
limited  size  of  the  plate.   With  only  three  equations  large 
errors  can  be  expected.   The  errors  can  be  greatly  reduced 
if  more  than  three  equations  are  generated  for  each  object 
point,  and  a  linearly  overdetermined  least  squares  method 
is  used  to  solve  them.   Subroutine  LLSQ  in  the  IBM  library 
was  used.   Nine  equations  were  used  to  determine  the  three 
unknown  displacements. 

The  accuracy  of  the  measurements  of  the  physical  set- 
up was  within  0.1  inch  and  the  accuracy  of  the  fringe  count 
was  estimated  to  be  +  0.1  fringe  number.   See  Table  II  for 
a  sensitivity  test. 

A  further  refinement  of  Equation  (7)  is  to  divide 
through  by  L.   The  three  unknowns  then  become  nondimensional 
displacements . 

L  '      L  '      L 

Thus,  even  if  the  wave  length  of  the  laser  is  unknown,  a 
relative  displacement  field  can  be  obtained. 

The  displacement  vector  obtained  from  LLSQ  may  or 
may  not  have  to  be  multiplied  by  (-1) .   The  sign  of  the 
displacement  vector  cannot  be  determined  from  holographic 
measurements  alone  because  the  reconstructed  image  and  the 
direction  of  fringe  movement  remains  unaffected  even  if  the 
direction  of  the  displacement  vector  is  reversed  by  revers- 
ing the  sequence  of  the  two  exposures.   A  priori  knowledge 
about  the  direction  of  at  least  one  displacement  component 
is  therefore  necessary  in  order  to  assign  proper  signs  to 

the  other  two  components. 
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If  the  fringe  distortion  about  a  point  is  severe, 
there  may  be  larger  than  0.1  fringe  count  errors.   A  tech- 
nique proposed  by  Dhir  to  detect  this  involves  plotting 
fringe  difference,  (nk~nk+l^  =  An,  along  a  given  cross 
section  of  the  object.   A  deformed  object  should  shov/  con- 
tinuous and  smooth  fringe  differences.   If  the  plots  of  A n 
are  not  smooth,  a  smooth  curve  is  passed  through  the  data. 
In  this  way  the  ZAn  for  a  distorted  fringe  area  can  be 
obtained  with  more  accuracy. 

2 .   The  Unfilled  Sample  Experiment 

The  sample  was  extended  to  0.135  inches  and  allowed 
to  relax  for  more  than  one-half  hour.   The  first  exposure 
was  then  made  and  the  wheel  was  turned  forty  degrees  to  make 
the  second  exposure  (Figure  12) .   It  was  found  by  experiment 
that  for  these  samples  a  forty-degree  turn  of  the  wheel  pro- 
duced the  best  fringe  density.   A  larger  turn  resulted  in 
the  fringes  being  too  close  for  accurate  counting.   A  smaller 
turn  produced  fringes  that  were  too  broad  and  increased  the 
difficulty  of  estimating  the  fringe  number  at  a  point. 

Nine  viewing  ports  were  drawn  on  the  holographic 
plate  with  india  ink  and  numbered  one  through  nine.   (See 
Figure  11) . 

The  location  of  the  zero  order  fringe  was  unknown, 
so  two  observations  of  a  single  grid  point  had  to  be  com- 
bined to  produce  one  equation.   Each  grid  point  was  observed 
nine  times  to  ensure  enough  over-determination  for  the 
least  squares  computer  solution.   The  nine  viewing  ports 
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were  used  in  the  following  combinations  for  each  point: 
1-7,  1-8,  1-9,  2-7,  2-8,  2-9,  3-7,  3-8,  3-9  or  1-3,  1-6, 
1-9,  4-3,  4-6,  4-9,  7-3,  7-6,  7-9.   If  certain  combinations 
were  unreadable,  variations  on  the  above  were  used. 

For  example,  when  determining  the  fringe  difference 
for  set  1-7,  the  object  grid  point  would  be  sighted  through 
port  number  one  and  assigned  an  arbitrary  fringe  number. 
If  the  grid  point  viewed  through  port  one  was  half  way  be- 
tween two  dark  fringes,  it  would  be  assigned  a  fringe  number 
of  0.5.   The  line  of  sight  was  then  slowly  changed  to  view- 
port seven  while  keeping  count  of  how  many  fringes  passed 
over  the  grid  point.   If  the  grid  point  when  viewed  through 
port  seven  was  between  fringes  three  and  four,  then  the 
number  3.5  would  be  assigned.   For  purposes  of  the  LLSQ 
solution,  the  fringe  difference  would  be  3.0.   The  actual 
tenth  of  a  fringe  number  assigned  for  each  point  was  deter- 
mined by  visual  interpolation  between  centers  of  dark 
fringes.   Most  of  the  fringes  could  be  readily  estimated  to 
plus  or  minus  one-tenth  of  a  fringe.   For  crowded  or  severely 
distorted  fringes  this  was  not  always  possible.   This  is  why 
the  fringe-difference  curve  smoothing  was  used. 

In  addition  to  the  59  grid  points  used  there  were 
seven  defined  points  on  each  crosshead  to  provide  a  refer- 
ence for  the  determination  of  rigid  body  displacements  and 
the  true  extension  of  the  sample. 

The  smoothed  data  were  used  to  solve  for  the  total 
U,V,  and  W  of  each  grid  point.   The  U,  V,  and  W  obtained 
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consisted  of  two  parts,  the  displacement  relative  to  the 

axis  system  and  the  displacement  due  to  rigid  body  affects. 

Thus: 

U  =  Ud  +  Urb 
V  =  Vd  +  Vrb 
W  =  Wd  +  Wrb 

The  subscript  d  refers  to  displacements  relative  to 
the  axes.   This  is  the  quantity  desired  for  comparison  with 
the  computer  solution.   The  subscript  rb  refers  to  the  rigid 
body  displacement.   The  rigid  body  displacements  were  deter- 
mined and  subtracted  from  the  total  displacements  to  yield 
the  relative  displacements. 

The  rigid  body  displacements  of  the  sample  were 
determined  from  the  displacement  of  the  X  and  Y  axes.   It 
was  assumed  that  the  X  and  Y  axes  remained  straight  even 
under  stress  due  to  the  geometric  symmetry.   Therefore,  any 
V  displacement  along  the  X-axis  was  due  completely  to  rigid 
body  displacement.   Any  U  displacement  along  the  Y-axis  was 
due  completely  to  rigid  body  displacement.   That  is,  U  =  Urb 
along  the  Y-axis  and  V  =  Vrb  along  the  X-axis. 

As  stated  previously,  a  priori  knowledge  of  one 
displacement  component  must  be  known  before  the  proper  sign 
can  be  given  to  the  other  two  components.   It  was  known  that 
the  V  displacements  of  the  nodal  points  on  the  Y-axis  were 
positive  since  that  was  the  sign  of  the  upper  crosshead 
extension.   Therefore  the  sign  of  the  U  and  W  displacements 
was  also  determined  for  the  Y-axis.   A  linear  least  squares 
curve  was  fitted  through  the  U  displacements,  making  Urb  a 
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function  only  of  Y.  It  was  then  assumed  that  the  Urb  dis- 
placement of  any  point  on  the  body  was  a  function  of  its  Y 
coordinate  only. 

By  letting  Y  equal  zero,  the  Urb  for  points  on  the 
X-axis  was  determined.   Although  Ud  along  the  X-axis  was 
unknown  its  sign  was  known  to  be  negative  from  the  Poisson 
effect  of  a  positively  strained  body.   Therefore, 

U  -  Urb  ^  0 
and  in  addition  Ud  must  become  more  negative  as  the  X 
distance  from  the  origin  increases. 

As  an  example,  the  following  is  a  sample  of  the  data 
obtained  for  two  nodal  points  on  the  X-axis.   (Figure  6) 

Point         U  V  W 

3        8.2512      3.8745      30.7820 
6       11.2715      4.0733      28.3705 

All  the  values  are  positive  but  it  was  not  known  whether 
they  should  be.   It  was  determined  using  the  method  described 
that  Urb  for  the  X-axis  was  -2.5057.   Applying  the  formula 
to  point  three, 

±8.2512  -  (-2.5057)  ^  0 
shows  that  the  U  for  point  three  must  be  negative  and  there- 
fore U  and  W  are  also  negative.   Then  Ud3  =  -5.7455  and 
Ud6  =  -8.7658  and  the  second  requirement  is  met.   With  the 
Vrb  for  the  X-axis  known  to  be  negative,  these  points  were 
also  fitted  with  a  linear  least  squares  curve  and  it  was 
assumed  that  the  Vrb  of  any  point  was  a  function  of  its  X 
coordinate  only. 
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The  determination  of  Wrb  is  very  difficult  and  can 
only  be  approximated  at  best.   There  are  no  sections  of 
the  sample  that  can  be  assumed  to  remain  straight  for  dis- 
placements in  the  Z  direction.   In  addition,  the  maximum 
Wd  displacement  was  predicted  from  the  computer  solution  to 
be  only  on  the  order  of  two  wave  lengths  with  most  points 
being  on  the  order  of  one  or  less,  which  borders  on  the 
expected  accuracy  limits  of  the  fringe  count  displacement 
prediction  method.   For  this  reason  W  was  not  considered 
in  the  comparisons  of  displacements. 

As  a  point  of  interest,  it  was  noted  during  the 
counting  of  fringes  that  the  maximum  fringe  shift  would 
occur  if  the  observation  viewports  used  were  in  a  line 
parallel  to  the  maximum  displacement  component. 

For  the  unfilled  material  experiment  the  additional 
extension  of  the  upper  end  was  found  to  be  24.7041  wave 
lengths  of  light.   Forty  degrees  of  handwheel  turn  was 
expected  to  cause  22.3  wave  lengths  of  additional  extension. 
3.   The  Filled  Sample  Experiment 

The  methods  used  to  make  the  interf erograms  and 
obtain  the  data  for  the  filled  sample  were  identical  to  those 
of  the  unfilled  sample.   The  only  change  was  in  the  amount  of 
extension  given  the  filled  sample  prior  to  the  making  of  the 
interf erogram.   It  was  found  by  experiment  that  the  samples 
would  fail  by  cracking  when  the  extension  was  about  0.225 
inches.   The  sample  was  extended  to  0.2  0  inches  to  ensure 
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that  dewetting  would  take  place.   Figure  13  is  a  photograph 
of  the  fringe  pattern  obtained. 

For  the  filled  material  experiment  the  additional 
extension  of  the  upper  edge  was  found  to  be  18.9060  wave 
lengths  of  light. 


24 


IV.   OTHER  METHODS 

A.  SIMPLE  EXTENSION 

The  simplest  method  attempted  was  to  apply  a  small  ex- 
tension between  exposures  and  visually  search  for  anomalies 
in  the  fringe  pattern.   No  such  pattern  was  detected  which 
could  be  attributed  to  a  dewetted  zone.   Figure  13  is  a 
typical  example  of  a  fringe  pattern  obtained.   Unwanted 
rigid  body  motion  may  have  caused  fringe  patterns  which 
disguised  the  zone. 

B.  RESIDUAL  STRAIN 

The  residual  strain  method  makes  use  of  material  prop- 
erty changes.   The  object  was  photographed  at  one  state  of 
strain  and  then  the  crosshead  was  extended  another  0.005 
inch.   The  increased  strain  caused  more  of  the  sample  to 
dewet.   The  crosshead  was  then  returned  to  its  original 
position  and  the  second  exposure  made.   The  material  that 
changed  properties  between  exposures  will  not  be  able  to 
return  to  its  original  state  and  this  will  cause  fringes. 
Figures  14  and  15  are  examples  of  this  technique.   The 
resulting  fringe  patterns  vary  widely.   On  several  occa- 
sions, patterns  appeared  in  one  or  more  quadrants  that  may 
or  may  not  have  been  dewetted  zones.   Fringe  patterns  on  the 
crossheads  indicated  that  they  were  not  returned  to  exactly 
the  original  extension  and/or  that  rigid  body  motion  had 
taken  place  between  exposures. 
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C.   TIME  AVERAGE 

Samples  of  the  filled  material  were  acoustically  excited 
and  time  average  holograms  were  made.   The  range  of  ex- 
citation frequencies  used  was  from  60  to  1000  cycles  per 
second.   The  holograms  obtained  showed  no  indication  of  a 
dewetted  area.   (Figure  16) 
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V.   RESULTS  AND  CONCLUSIONS 

A.   THE  UNFILLED  MATERIAL 

The  computer-predicted  displacements  and  the  measured 
displacements  for  the  unfilled  material  are  plotted  versus 
radial  distance  for  each  ten-degree  radial  of  the  finite 
element  grid.   (Figures  17  through  21) .   An  examination  of 
the  plots  indicates  that  the  linear  theory  does  not  accu- 
rately predict  the  behavior  of  the  sample.   The  measured 
U  and  V  displacements  are  higher  than  predicted  on  every 
radial . 

There  are  indications  that  the  deviation  from  linearity 
depends  upon  the  intensity  of  the  stress.   It  can  be  seen 
that  the  V  measured  points  match  very  well  with  the  shape 
of  the  V  predicted  curve  at  the  higher  radials.   Since  the 
displacement  at  any  point  is  the  integral  of  the  strains 
up  to  that  point,  it  would  indicate  that  the  material  at 
the  higher  radials  is  behaving  linearly  and  the  larger  than 
predicted  displacements  are  due  to  the  nonlinear  strains  on 
the  lower  radials. 

Another  unexpected  result  is  that  the  measured  U  for 
the  ten-degree  radial  is  larger  than  that  of  the  zero-degree 
radial.   The  U  for  the  20-degree  radial  is  smaller  than  the 
U  for  the  10-degree  radial  but  still  larger  than  the  zero. 
It  is  not  until  the  30-degree  radial  that  the  measured  U  is 
smaller  than  that  of  the  zero-degree  radial.   The  predicted 
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U  is  greatest  for  the  zero-degree  radial  and  decreases  for 
each  increase  in  radial. 

B.  THE  FILLED  MATERIAL 

For  the  filled  material  the  linear  predictions,  the 
dewetted  predictions  and  the  measured  displacements  are 
plotted  versus  radial  distance  for  each  radial.   (Figures 
22  through  26) . 

The  linear  V  displacements  predictions  were  in  general 
better  than  the  dewetted  predictions.   This  is  due  in  part 
to  the  early  prediction  of  failure  from  the  dewetted  pro- 
gram.  The  failure  prediction  was  discounted  and  the  result 
was  that  the  predicted  dewetted  zone  grew  to  cover  most  of 
the  right-hand  side  of  the  grid.   The  lowered  elastic  moduli 
of  the  dewetted  lower  radials  allowed  too  much  strain  to 
occur  in  the  elements  along  the  X-axis,  with  the  resultant 
larger  displacements  at  the  higher  radials. 

Both  computer  methods  underestimate  the  U  displacements 
and  they  agree  with  each  other  very  well  for  the  30-degree 
radial  and  higher. 

The  measured  values  of  V  fell  between  the  linear  and 
dewetted  computer  solutions  but  were  closer  to  the  linear 
solution.   This  indicates  that  the  material  did  dewet  but  to 
a  much  smaller  extent  than  predicted  by  theory. 

C.  FILLED  DISPLACEMENTS  COMPARED  WITH  UNFILLED  DISPLACEMENTS 
The  measured  displacements  for  the  filled  and  unfilled 

material  are  plotted  versus  radial  distance  for  each  ten 
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degree  radial  (Figures  27  through  31).   The  unfilled  dis- 
placements were  converted  to  the  same  upper  edge  extension 
as  the  filled  material  by  multiplying  by 

18.9060  _  n  ^„0fi 
24.7041  '  °-765298 

The  unfilled  displacements  are  greater  than  the  filled 
displacements  for  every  radial  but  in  general  the  shapes 
of  the  displacement  curves  are  similar.   This  appears  to  be 
due  to  nonlinear  straining  of  the  unfilled  material  in  the 
lower  radials. 

D.  THE  HOLOGRAPHIC  DISPLACEMENT  DETERMINATION  METHOD 
The  method  for  obtaining  general  displacements  from  a 

holographic  interf erogram  proved  to  be  simple  to  understand 
and  easy  to  use.  However,  its  use  is  not  recommended  for 
general  use  because  of  the  large  amount  of  time  that  must 
be  spent  to  obtain  data.  It  is  estimated  that  a  total  of 
twelve  hours  were  spent  just  to  determine  the  nine  fringe 
differences  for  each  of  the  73  points  used  in  the  analysis. 

E.  THE  TENSILE  MACHINE 

The  tensile  machine,  as  designed,  was  not  rigid  enough. 
The  unwanted  movement  of  the  machine  induced  fringe  patterns 
onto  the  sample  that  disguised  the  true  patterns  and  made 
simple  visual  inspection  of  the  patterns  impossible.   Any 
continuation  of  this  work  should  be  done  with  a  much  more 
rigid  tensile  testing  machine. 
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F.  EQUIVALENT  STRESS 

The  concept  of  equivalent  stress  was  derived  for  use 
in  predicting  plasticity  in  metals.   It  is  used  to  relate 
a  three-dimensional  stress  state  to  a  one-dimensional  one 
such  as  a  stress-strain  diagram,  however  it  is  only  a 
postulate,  not  a  proven  fact.   The  effects  of  plasticity 
and  dewetting,  as  far  as  the  stress-strain  diagrams  are 
concerned,  are  the  same  so  it  seemed  a  logical  extension 
of  the  equivalent  stress  idea  to  apply  it  to  dewetting  of 
propellants.   The  results  indicate  that  this  is  not  the 
case  and  that  equivalent  stress  may  not  be  applied  to 
dewetting  with  accurate  results. 

G.  THE  FINITE  ELEMENT  PROGRAM 

From  the  inaccuracy  of  the  computer  displacement  pre- 
dictions and  the  early  prediction  of  propellant  sample 
failure  it  is  evident  that  a  finite  element  material 
property  change  criterion  based  on  equivalent  stress  alone 
will  not  give  satisfactory  results.   The  effects  of  visco- 
elastic  stress  relief  will  have  to  be  incorporated  into  the 
computer  dewetting  criterion  so  that  stresses  will  build 
more  slowly  without  over-predicting  the  resultant 
displacements. 

The  linear  finite  element  program  more  closely  predicted 
the  V  displacements  of  the  propellant  but  it  also  predicted 
a  stress  of  over  200  psi  at  node  one  for  a  0.20  inch  extension, 
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which  is  over  twice  the  failure  stress  obtained  from  the 
tensile  test.   The  sample  did  not  fail  at  0.20  inches  of 
extension  which  emphasizes  the  importance  of  viscoelastic 
stress  relief  incorporation  into  the  finite  element  model. 

H.   CONCLUDING  REMARKS 

The  comparison  of  the  computer  predicted  displacements 
to  the  measured  data  was  unsuccessful  due  to  the  inconsist- 
ency of  the  results.   Modifications  to  the  computer  tech- 
nique as  discussed  in  section  G  should  produce  better 
results.   In  addition,  a  more  rigid  tensile  testing  machine 
would  eliminate  errors  resulting  from  rigid  body  motion  and 
unsymmetric  crosshead  extension. 

The  dewetted  zone  was  not  detected  visually.   The  most 
promising  technique  appears  to  be  the  residual  strain  method 
and  further  investigation  of  the  technique  is  recommended. 
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TABLE  I 
TABULATION  OF  MATERIAL  PROPERTIES 

MATERIAL  TYPE          POISSON'S  RATIO    ELASTIC  MODULUS 

1  0.485  1091.0 

2  0.483  905.7 

3  0.481  754.7 

4  0.480  294.3 


32 


TABLE  II 

TEST  of  SENSITIVITY  of  LEAST  SQUARES  FRINGE  DIFFERENCE 

METHOD 


Data  used: 

Object  Number 

First 

View 

Port 

Second 

Po 

rt 

An 

1 

1 

3 

3.7 

1 

1 

6 

3.8 

1 

1 

9 

3.1 

1 

4 

3 

4.0 

1 

4 

6 

,   4.1 

1 

4 

9 

3.4 

1 

7 

3 

4.3 

1 

7 

6 

4.4 

1 

7 

9 

3.7 

Test  results: 


Numbe 

r  of 

Solution 

Equat 

ions 
4 

used 

U 

V 

W 

-7. 

633 

-3. 

,172 

-19. 

69 

5 

-7. 

637 

-3. 

,437 

-21. 

83 

6 

-7. 

875 

-3. 

,020 

-17. 

,27 

7 

-7. 

894 

-2, 

,489 

-14. 

,38 

8 

-7. 

896 

-2. 

,542 

-15. 

,07 

9 

-7. 

829 

-2. 

,544 

-15. 

,07 

9 

-7. 

799 

-2. 

,550 

-15. 

19 

note 

(1) 

9 

-8. 

053 

-2. 

,544 

-15. 

,08 

note 

(2) 

9 

-7. 

605 

-2. 

,544 

-15. 

,06 

note 

(3) 

Note  (1)   Odd  numbered  data  cards  had  0.1  subtracted 

from  /\n  and  even  numbered  cards  had  0.1  added. 

Note  (2)   All  nine  data  cards  had  0.1  added  to  An. 

Note  (3)   All  nine  data  cards  had  0.1  subtracted  from  An 
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Figure  3 
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Finite  Element  Grid 
Figure  5 
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Figure  6 
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Figure  9 
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Time  Average  Hologram  of  Filled  Sample  Excited  at  120  CPS 

Figure  16 
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Figure  17 
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